Dynamics and Transport in Random Antiferromagnetic Spin Chains. 
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We present the first results on the low-frequency dynamical and transport properties of random 
antiferromagnetic spin chains at low temperature (T). We obtain the momentum and frequency 
dependent dynamic structure factor in the Random Singlet (RS) phases of both spin-1/2 and spin-1 
chains, as well as in the Random Dimer phase of spin-1/2 chains. We also show that the RS phases 
are unusual 'spin-metals' with divergent low-frequency conductivity at T = 0, and follow the spin 
conductivity through 'metal-insulator' transitions tuned by the strength of dimerization or Ising 
anisotropy in the spin-1/2 case, and by the strength of disorder in the spin-1 case. 
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Low-dimensional quantum spin systems in the pres- 
ence of disorder are a fascinating laboratory for the study 
of the interplay between quantum effects, strong corre- 
lations and disorder. This is in part because they are 
often very sensitive to disorder, and display dramatic 
effects due to small amounts of disorder (especially in 
one dimension). The question of disorder effects is thus 
of considerable experimental importance in such cases. 
Moreover, a great deal is now known about the effects 
of strong correlations and quantum fluctuations in pure 
one and two-dimensional spin systems |Q, allowing one 
to focus on the new features introduced by disorder. 

Random antiferromagnetic spin-1/2 chain compounds 
are interesting, experimentally realizable examples [^]|| 
of such systems. Theoretical work has led to the pre- 
diction that even a small amount of disorder can cause 
these systems to have an extremely unusual 'Random 
Singlet' (RS) ground state for a range of parameter val- 
ues In this RS state, the interplay of disorder and 
quantum mechanics locks each spin into a singlet pair 
with another spin; the two spins in a given singlet pair 
can have arbitrarily large spatial separation. A similar 
RS state is also predicted in the spin-1 Hcisenberg an- 
tiferromagnetic chain for sufficiently strong randomness 
in the bonds HJ7|]. The theory also yields unusual re- 
sults for the low-temperature thermodynamics and some 
ground-state correlators in the RS states. 

In this Letter, we report on the first theoretical study 
of the dynamical and transport properties of these sys- 
tems; the primary motivation is of course to understand 
the dynamics and transport when quantum interference, 
strong correlations and disorder are all simultaneously 
important. We focus here on the momentum and fre- 
quency dependent dynamical structure factor S(k, uj) and 
the dynamical spin conductivity. The former is directly 
probed by inelastic neutron scattering (INS) experiments 
on these compounds, while the latter is of considerable 
theoretical interest as it contains information about the 
nature of any localization phenomena, or lack thereof, in 
these strongly correlated disordered systems. 



In particular, we calculate the low-temperature S(k, uS) 
at low frequencies in the RS states of spin-1/2 and spin- 
1 chains as well as in the 'Random Dimer' (RD) phase 
of spin-1/2 chains that is predicted ||] in the presence 
of weak enforced bond-dimerization and strong disorder. 
Our results lead us to expect some very unusual features 
in the INS cross-section as a function of wavevector trans- 
fer k (for fixed, small energy transfer oj > T), especially 
in the RD phase. For the conductivity, we calculate the 
real part of the low-frequency dynamical spin conduc- 
tivity a'(ui) at low temperatures and show that the RS 
phases of both spin-1/2 and spin-1 chains are actually un- 
usual 'spin-metals' with a <j'(u)) that diverges for low lu at 
T = 0. We also follow <j'(lu) across novel 'metal-insulator' 
transitions to insulating states; these are accessed in the 
spin-1/2 case by turning on dimerization (leading to a 
RD phase) or by increasing the Ising anisotropy above a 
threshold, and in the spin-1 Heisenberg case by reducing 
the disorder below a critical value. Thus, we are able to 
obtain a wealth of reliable information about the unusual 
dynamics and transport in these strongly correlated ran- 
dom quantum systems. We emphasize that all our results 
are exact in the low- frequency limit at T = 0, and con- 
tinue to be valid for T ^ so long as T < u) (in some 
cases, we are also able to access frequencies lu < T pof). 

In the spin-1/2 case, the specific problem we consider 
is the random XXZ Hamiltonian, which describes the low 
energy (magnetic) dynamics of insulating antiferromag- 
netic spin-1/2 chain compounds with chemical dis- 
order that affects the bond strengths: 
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where Sj are spin-1/2 operators at lattice sites j sepa- 
rated by spacing a, and both jj- and JJ are random 
positive exchange energies with a joint probability distri- 
bution Pq(J ± , J z ). Detailed information on the nature of 
the excitations in such systems is encoded in the dynam- 
ical structure factor S al3 (k, lu) (where a (3 = H — or a (3 = 
zz) with the T = spectral representation 
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where s ± = s 1 ± is y and {|m)} denote the exact eigen- 
states of the system with energies E m (|0) is the ground 
state, with Eq = 0). 

As mentioned earlier, we characterize transport in 
terms of the dynamical conductivity <t(uj). The real 
part <j'{lo) of <j(uj) is defined by the relation P(u) — 
a'(uj)\VH\ 2 (iv), where P{oj) is the power absorbed per 
unit volume by the system when a uniform magnetic field 
gradient VH(lo) (where the field H always points in the 
z direction) oscillating at frequency u is applied along 
the length of the chain. From standard linear response 
theory, we have the following Kubo formula for <t'(u>) at 
T = 0: 
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s^sJ +1 )/2 is the current op- 
erator on link j that transfers one unit of the z com- 
ponent of the spin from one site to the next, and the 
frequency uj is taken positive for notational convenience. 
Note that both S a/3 (k,u>) and cr'(w), as defined here, are 
self-averaging in the thermodynamic limit. 

Randomness in the bonds is a relevant perturbation |j| 
for the pure XXZ chain when < J z /J ± < 1; any 
amount of disorder is thus expected to drive the system to 
strong disorder. In this regime, the system can be treated 
by a strong-randomness renormalization group (RG) that 
proceeds as follows : We look for the bond with the 

largest J 1 - in the chain, say J^ 3 between spins 2 and 3 
- this sets the energy cutoff = maxjjj 1 }. We first 
solve the corresponding two spin problem (the neigh- 
bouring bonds are introduced later as perturbations). So 
long as the J z are not large compared to the J , the 
ground state of the two-spin problem will always be a 
singlet separated by a large gap from the triplet excited 
states. We can then trade our original Hamiltonian in for 
another Hamiltonian (determined perturbatively in the 
ratio of the neighbouring bonds to the strongest bond) 
which acts on a truncated Hilbert space with the two 
sites connected by the 'strong' bond removed. To lead- 
ing order, this procedure renormalizes the Hamiltonian 

whites = Ej 3 =i JH^$+i +a i 8 i+i) +J ? 8 3 a l+i to7 ^ 14 = 
Ji(sfjt + 44) + k44 with Ji = JirJi/iJi + Ji) 

and Jf = JfJ§/2J^- ||. Note that the new bond has 
length li = h+h + h- This procedure, if it remains valid 
upon iteration, thus ultimately leads to a ground state of 
singlet pairs, with pairs formed over long distances held 
together by correspondingly weak bonds; this is the RS 
state alluded to earlier. 

A complete understanding of the possible states thus 
requires an analysis of the effects of iterating the ba- 
sic RG procedure. Such an analysis was performed in 



Ref H leading to the following conclusions: If the J z 
dominate over the J -1- , this procedure rapidly becomes 
invalid and the ground state actually has Ising Antifer- 
romagnetic (IAF) order. In all other cases, the ground 
state is a Random Singlet state. The low-energy (with 
energy cutoff <C O c , where O c is the microscopic cut- 
off) effective theory in these cases is written in terms of 
the n(£l)L 'surviving' spin variables (the bond-length I 
between successive surviving sites is now a random quan- 
tity). When the J- 1 dominate, this effective Hamiltonian 
has all J z = 0, and J and I drawn from a universal 
joint probability distribution V(J ± ,l\fl) characteristic of 
the 'XX Random Singlet fixed point' (XX-RS) of the 
RG. Between the IAF phase and this XX-RS phase lie 
two kinds of critical points. If the initial problem has 
full Heisenberg symmetry (J z = J- 1 for each bond), the 
low-energy effective Hamiltonian preserves this symme- 
try and has bonds strengths and lengths drawn from the 
same probability distribution V(J,l\Cl). In the RG lan- 
guage, the Heisenberg system is controlled by the XXX- 
RS fixed point. Finally, in this language, the generic crit- 
ical state is controlled by the XXZC-RS fixed point — the 
low energy effective theory has bonds and lengths drawn 
from a fixed point distribution V\{J , J z ,l\£l) with the 
property / dJ z Vi = T>{J^,l\VL). The probability distri- 
butions V and V\ become infinitely broad as f2 — * 0; 
this implies that the RG becomes asymptotically exact 
at low-energies, and in particular predicts the ground 
state properties and low-temperature thermodynamics 
correctly. Note that in all of the above, we have sup- 
pressed non-universal scale factors multiplying the argu- 
ments of the functions n and V; these scaling functions 
however remain the same for all systems that flow to any 
of the RS fixed points. 

To use the foregoing for the calculation of dynamical 
or transport properties, we need to explicitly keep track 
of the renormalization of the operators that enter spec- 
tral sums such as (0) and (S). We illustrate our general 
approach jlO) by first calculating the T = dynamical 
conductivity at asymptotically low frequencies. The first 
step is to work out the rules that govern the renormal- 
ization of the current operators. Assume, once again, 
that J^j is the strongest bond. We wish to work out 
perturbatively the renormalized operators f 1/2/3 that we 
trade in ri/2/3 for, when we freeze spins 2 and 3 in their 
singlet ground state (the other current operators to the 
left and right of this segment of our system are left un- 
changed to leading order by the renormalization). Now, 
note that these other operators have an overall scale fac- 
tor in them that is nothing but the corresponding J . 
In order to be consistent, we clearly need to work out 
T1/2/3 correct to 0(Jj L ) by adding the effects of virtual 
fluctuations to the projection of 71/2/3 m to the singlet 
subspace 1 1 . An explicit calculation gives the sim- 
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three operators renormalize to the same operator, which 
we will denote henceforth by f\ for consistency of nota- 
tion. 

As we carry out the RG and reduce the energy cut- 
off, the above result implies that Ylj T j renormalizes to 

yijljfj, where j now labels the uqL sites of the renormal- 

ized lattice at cutoff scale il, and the lj are the lengths 
of the renormalized bonds in this problem. With this 
in hand, we run the RG until the cutoff is reduced to 
Q u = lo (note that the operator fj linking two sites con- 
nected by a strong bond {Jj~, Jj} promotes the corre- 
sponding pair of spins from their singlet ground state to 
the triplet state | to) with m z = 0, which is separated 
from the ground state by precisely Jj 1 ) and rewrite @ 
as 



loL 



(4) 



where the tildes are a reminder of the fact that this spec- 
tral sum now refers to the new Hamiltonian with cutoff lo 
which has only n u L sites (the couplings and bond lengths 
in this problem are of course drawn from the probability 
distribution characteristic of the fixed point to which the 
system flows in the low energy limit). 

We now have to calculate the spectral sum (^) in this 
new problem. The following crucial observation allows us 
to do this: At the next step of the RG, one would have 
looked for all the bonds in this renormalized problem 
that have J 1 - in the range (fi^ , f2 w — dil) and formed sin- 
glets out of the corresponding pairs of spins. The states 
with E m = lo that give the dominant contribution to the 
spectral sum (Q) correspond precisely Jl2| to promoting 
any one of these pairs to the triplet state |io) under the 
action of the current operator living on the correspond- 
ing bond. The matrix element for this transition is just 
Ilo/2, where I is the length of the bond in question. In 
the thermodynamic limit, we thus have 

^(oA^I^k f dldJ^u 2 l 2 V{J^,l\n u )8(u - J x ) 
w J 

~ ln(n c /w) ; (5) 

the last line is the leading behaviour for lo <C Q c obtained 
by using the results of Ref J(| for nn and V . Notice 
that this analysis holds equally well at all three RS fixed 
points. 

Thus cr'(w) diverges logarithmically for small enough lo 
in the unusual 'spin-metal' phase controlled by the XX 
fixed point [ p"3"| as well as at the critical point separat- 
ing this phase from the insulating phase with IAF or- 
der in the ground state. Note that this 'metal-insulator' 
transition has the curious feature that the quantum crit- 
ical point separating the conducting phase from the IAF 
insulator has the same T = transport properties as 
the conducting phase. On the insulating side, o'{lo) is 



suppressed below a pseudo-gap energy E g . The domi- 
nant contributions for to <C E g come from Griffiths ef- 
fects in which rare fluctuations in the couplings of the 
Hamiltonian allow a long finite segment of the system 
to be 'locally' in the 'metallic' phase. A calculation of 
this contribution gives [jl0| a low-frequency conductiv- 
ity o~' {lo) ~ Lo a ln 2 (f2 c /(jj), where a > is a continu- 
ously varying, non-universal exponent that vanishes at 
the transition. 

The dynamical structure factor can be calculated in a 
similar way. Consider first S zz (k,Lo). One begins with 
the spectral sum (|^). The leading order 'operator renor- 
malizations' needed in this case are particularly simple — 
each spin operator remains unchanged so long as it is 
not part of a singlet and renormalizes to zero upon being 
locked into a singlet state As before, we run the 

RG until the cutoff = lo and do the spectral sum with 
the renormalized operators in the new problem. This 
renormalized sum may be evaluated by again recogniz- 
ing that it is dominated by excitations to the triplet 
state \to) of pairs of spins connected by the renormalized 
bond J 1 - = lo. The corresponding matrix element is sim- 



ply (i 
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)/2, where I is the length of the strong bond 



(note that If a is an odd integer). This allows us to write 
for lo <C tt c and k = tt / a + q 

S zz {k,Lo) ~7i n „ / dldJ x \l + e iql \ 2 V(J X ,l\u)6{u;- J x ) 



at all three RS fixed points. The calculation of (k, lo) 
is slightly more involved as the gap to the relevant triplet 
excited state (with m z = 1) of a pair of spins connected 
by a strong bond (J ± , J z ) is now (J 1 - + J z )/2. However, 
a careful analysis (!(]] gives the same result as above for 

LO <C Q c . 

Let us focus here on the regime \q\ = \k — ir/a\ <C a -1 
(in addition to lo <C fi c ). In this regime, the integral can 
be evaluated using the results of Ref [§ for V . This gives 
the following rather unusual universal scaling form at the 
RS fixed points (i.e so long as the ground state does not 
have IAF order) 
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where a(3 = H or zz, A and are non- universal 

scale factors, and the fully universal function $(a;) can 
be written as 



$(x) = 1 



cc(cos(a;) sinh(x) + sin(x) cosh(a;)) 
cos 2 (x) sinh 2 (x) + sin 2 (a;) cosh 2 (x) 



(7) 



We now briefly summarize the effects of weak (com- 
pared to the disorder) enforced bond dimerization con- 
trolled by a dimensionless small parameter |<5| that sets 
the difference between the probability distributions of 
even and odd bonds [Q. The RG flows in the vicin- 
ity of the XX and XXX RS states are known ||[|; the 
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low energy properties are controlled by lines of strong- 
disorder Random Dimer || fixed points ending in the 
XX and XXX fixed points. These RD fixed points again 
describe ground states that consist of singlet pairs; how- 
ever, now the singlet bonds preferentially start from an 
even (odd) site and end at an odd (even) site for positive 
(negative) 8. Our approach readily allows us to follow the 
full crossover of the low-frequency dynamical conductiv- 
ity and structure factor from the XX and XXX-RS states 
to the corresponding RD phases (only the non-universal 
prefactors differ in the two cases). Here, we focus on re- 
sults deep in the RD phases, i.e for frequencies such that 
T u = \5\ \n(Q c /uj) ^> 1. For the conductivity, we obtain 
c 1 {lo) ~ ln(ri c /tj)r [i; e _Cor " , which can be wriiten more 
explicitly as a' (w) - \S\uj C0 ^\n 2 (fl c /u;), where Co is a 
non-universal constant. The RD phases are thus seen to 
be gapless insulators. The dynamic structure factor in 
the vicinity of k = it j a can be written as 

S °0(k = - + q,oj) = ^(1 + cos( Cl r>-)e- C2f ^ 2 ), (8) 
a to 1 a 

where a = cq\S\, ci, C2 and C are non- universal constants, 
and q = qa/S 2 is assumed <C 1. This result has striking 
oscillatory structure which is best understood |l^] as a 
novel signature of the sharply defined geometry of the 
rare Griffiths regions that contribute to the scattering 
at a given low frequency — more precisely, the average 
length of the relevant Griffiths regions is of order r w a/|<5| 2 
while the RMS fluctuations in the lengths are only of 

order v^WI^I 2 - 

In the spin-1 case, the specific Hamiltonian we consider 

is 

n = J2JjS r s J+1 , (9) 

3 

where Sj are spin-1 operators and the Jj are random, 
positive bond-strengths; the corresponding distribution 
of ln(Jj ) is characterized by a width W. A renormaliza- 
tion group analysis reveals that the system flows to the 
analog of the XXX RS point described earlier only when 
W exceeds a critical value W c . In this case, our previous 
results for the dynamic structure factor and the dynam- 
ical conductivity continue to apply. As W is decreased, 
the system undergoes a quantum phase transition to the 
so-called Gapless Haldane (GH) phase Q|; both the quan- 
tum critical point and the GH phase in the vicinity of 
it are however still controlled by strong-disorder fixed 
points 0. Our technique may also be used to calculate 
the dynamical conductivity at this critical point and in 
the GH phase, although the details differ considerably 
from the spin-1/2 case |k]]. At the critical point, we find 

a'(iu) ~ ln 2 {n c /co) , (10) 

which is a stronger divergence than in the strong-disorder 
RS phase. In the GH phase, the conductivity goes as 



o'{lu) ~ uj a ln 2 (f2 c /ijj) for log-frequencies ln(f2 c /u>) ^> 
(W c — W)~"/ 3 (y is the correlation length exponent 
known to equal 6/(-\/l3 — 1)). This phase is thus 
a gapless insulator, not unlike the RD phase of spin- 
1/2 chains. The non-universal exponent a is given as 

a ~ (W c - wy^. 

We close with some remarks regarding experiments. 
Our results for S{k, lj) near k — n/a are clearly of direct 
relevance to low-temperature INS experiments on such 
spin chain systems in the regime where the energy trans- 
fer lo satisfies T < lo <C O c . Of particular interest would 
be an experimental confirmation of the coherent oscil- 
latory structure predicted in the RD phase. Regarding 
transport, we hope that our results motivate experiments 
to probe the spin conductivity in these systems. 
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